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Abstract-  We  present  an  analysis  of  the  characteristics  of 
different  filters  for  the  compression  of  magnetic  resonance 
images.  Compression  rates  were  33:1  and  50:1.  We  compare  the 
performance  among  different  types  of  wavelets  presented  in  the 
literature  and  provide  quantitative  (percentage  of  energy 
retained,  peak  signal  to  noise  ratio)  and  qualitative  (analysis  by 
a  group  of  seven  experts)  data  to  support  our  conclusions. 
Different  types  of  coiflets,  symlets  and  biorthogonal  wavelets  are 
analyzed,  and  we  conclude  that  for  the  images  under  study  (T1 
weighed  images  in  three  planes),  the  best  results  are  provided  by 
the  biorthogonal  spline  (Daubechies)  wavelet  2,6.  Several 
explanations  for  these  results  are  mentioned. 

Keywords-  Image  compression,  wavelet  packets,  magnetic 
resonance  imaging 

I.  Introduction 

In  the  last  years  the  growth  of  digital  imaging  in  the  medical 
environment  has  outpaced  any  reasonable  expectations.  This 
growth  has  profited  from  the  developments  in  computer  and 
communications  technology  for  storage  and  transmission  of 
images  within  and  among  medical  establishments.  Even  with 
growing  bandwidth  and  diminishing  costs,  transmission  and 
storage  of  medical  images  quickly  runs  into  limitations.  One 
of  the  approaches  to  these  limits  is  to  compress  the  images  in 
order  to  reduce  both  storage  space  and  transmission  time. 
Although  several  image  compression  methods  such  as  IPEG 
have  been  standardized  for  some  time,  these  methods  have 
not  been  widely  accepted  in  the  clinical  environment.  This  is 
due  to  the  fact  that  in  order  to  achieve  high  compression 
rates,  inexact  or  "lossy"  algorithms  such  as  the  discrete 
cosine  transform  have  to  be  employed,  and  these  algorithms 
introduce  distortions  or  artifacts  that  may  affect  the 
diagnostic  quality  of  medical  images. 

The  theory  of  wavelet  analysis  has  proved  to  be  a  very 
important  development  in  the  search  of  more  efficient 
methods  of  image  compression.  In  particular,  the  fact  that 
wavelet  coefficients  indicate  changes  in  the  signal  under 
analysis  means  that  areas  of  small  changes  or  with  no 
changes  at  all  will  produce  small-valued  coefficients.  In  the 
case  of  image  compression,  generally  there  are  many  areas 
of  this  type  in  a  typical  image  so  image  compression  using 
this  scheme  should  be  quite  efficient. 

Different  properties  of  the  wavelets  have  been  studied 
under  the  point  of  view  of  efficient  compression.  In 
particular,  the  number  of  vanishing  moments  has  been  linked 
to  signal-to-noise  ratio  and  smoothness  v.s.  ringing  in 
compressed  image  reconstruction.  Other  properties  such  as 
symmetry  and  regularity  have  also  been  evaluated  in  this 
light  [1,2],  In  addition  to  this,  researchers  in  the  field  of 
vision  as  well  as  those  dedicated  to  wavelet  analysis  [3,4] 
have  proposed  a  link  between  the  properties  of  the  wavelets 
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response  of  the  human  visual  system  has  been  compared  to  a 
Gabor  function,  which  has  been  linked  to  some  wavelet  filter 
banks.  Because  of  this,  he  relationship  between  the 
amplitudes  of  the  primary  and  secondary  lobes  of  some 
wavelets  has  been  studied  in  the  analysis  of  the  performance 
of  different  wavelets  for  image  compression  [5].  Other 
authors  have  analyzed  quite  a  number  of  wavelet  filter  banks 
in  order  to  determine  those  with  the  best  characteristics  for 
compression.  In  this  respect,  some  of  the  most  cited  wavelets 
are  Coiflets,  which  have  the  highest  number  of  vanishing 
moments  for  a  given  support  width;  Symlets,  which  are 
nearly  symmetrical,  which  means  that  they  behave  like 
linear  phase  filters,  that  are  known  to  be  important  for  image 
compression,  and  biorthogonal  filters,  which  allow  for 
perfect  reconstruction  without  redundancy,  and  offer 
symmetry  and  compact  support  [6], 

The  purpose  of  this  paper  is  to  investigate  which  are  some 
of  the  parameters  that  are  directly  responsible  for  the  quality 
of  low  bit-rate  coded  images.  Several  filter  banks  that  have 
been  mentioned  in  the  literature  are  evaluated,  together  with 
parameters  such  as  filter  length,  symmetry  and  wavelet 
shape. 

II.  Methodology 

For  the  purpose  of  this  study,  four  types  of  image  were 
analyzed:  Three  T1 -weighed  magnetic  resonance  slices  in  the 
axial,  saggittal  and  coronal  planes,  and  a  common  Lena 
image.  Images  were  compressed  at  33:1  and  50:1  ratios  using 
the  wavelet  packet  decomposition  at  five  levels,  and  the 
Shannon  entropy  criterion  for  the  selection  of  the  best  tree. 
The  wavelet  coefficients  were  thresholded  in  order  to  achieve 
the  desired  compression  rate.  Coefficients  below  the 
threshold  were  set  to  zero.  Three  groups  of  wavelets  were 
employed  in  this  study;  Coiflets,  Symlets  and  Biorthogonal. 
A  preliminary  study  was  carried  out  in  order  to  select  the  best 
performing  wavelets  in  each  group.  The  wavelets  that  were 
retained  were  Symlets  6,  7  and  8;  Coiflets  2,3,  4  and  5; 
Biorthogonal  (Daubechies)  2,8,  3,9,  5,5,  and  6,8. 

Quantitative  analysis  consisted  in  the  measurement  of  the 
peak-  to-peak  ratio,  the  peak  signal  to  noise  ratio  and  the 
percentage  of  the  energy  retained  after  compression  for  both 
compression  rates,  as  well  as  the  determination  of  the 
symmetry  of  the  reconstruction  filter  T;. 

Qualitative  analysis  was  carried  out  by  a  group  of  seven 
experts  on  medical  images  from  the  Image  Processing 
Laboratory  at  UAM-Iztapalapa  and  the  National  Institute  of 
Neurology  and  Neurosurgery  at  Mexico  City.  Images  were 
displayed  on  high-quality  17"  flat  displays  set  at  1024  x  768 
pixel  resolution  at  an  85Hz  refresh  rate.  The  experts  graded 
the  images  on  a  10-point  scale,  ranging  from  0  (unusable 
images)  to  10  (images  deemed  equal  to  or  as  good  as 
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III.  Results 

Table  1  shows  the  data  for  the  quantitative  analysis  for  the 
four  images  and  all  the  filter  types.  The  second  column  shows 
the  values  for  the  percentage  of  the  energy  retained  after 
compression  (expressed  as  the  fraction  of  this  value  after 
99%).  The  middle  columns  show  the  peak  signal-to-noise 
ratio  for  the  two  compression  ratios  involved.  The  last  two 
columns  show  the  symmetry  and  the  peak-to-peak  ratios  for 
all  the  wavelet  types.  Qualitative  analysis  data  are  shown  on 
Figure  1.  The  mean  values  and  the  standard  deviations  are 
shown  for  all  filter  types. 
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Table  1.  Results  of  quantitative  analysis. 


Statistical  analysis  was  carried  out  first,  by  analysis  of 
variance,  where  the  null  hypothesis.  Ho  (all  the  means  are 
equivalent)  is  rejected  at  the  95%  confidence  level.  However, 
as  the  qualitative  grading  system  is  considered  to  be  a  "soft" 
variable,  the  condition  of  normality  cannot  be  assumed,  we 
carried  out  the  Friedman  test,  which  is  well-suited  for  an 
evaluation  where  different  judges  evaluate  several  classes  of 
objects.  In  this  case  Ho  (grades  are  not  correlated)  is  rejected, 
with  a  value  of  P  <0.00000001,  Q=48.09.  These  results  mean 
that  the  different  processes  can  be  distinguished  among  them. 
The  columns  have  different  values  (different  processes)  and 
the  rows  are  correlated  (there  is  agreement  among  grades).  In 
addition  to  this  study,  the  analysis  of  the  modes  of  the  grades 
together  with  the  average  values  of  the  grades  do  allow  us  to 
determine  which  is  the  best  type  of  filter  for  the  images  under 
study.  Figure  2  shows  the  different  reconstructions  for  the 
best  filters  of  each  of  the  different  major  types  of  wavelet. 
The  compression  ratio  in  this  case  is  50:1  in  order  to  make 
errors  more  noticeable. 

IV.  Discussion 

For  biorthogonal  wavelets,  b2.8  and  b3.9  the  percentage 
of  the  energy  recovered  from  the  50:1  compressed  images 
was  in  all  instances  over  99%,  and  this  occurred  only  for 
these  filters.  Other  filters  with  over  99%  of  the  energy 
recovered  were  b5.5  for  the  Lena  image  and  b5.5  and  b6.8  for 
the  axial  MR  images.  It  is  worth  noting  that  for  the  cases  of 
filters  b5.5  and  b6.8,  where  decomposition  and  reconstruction 
filters  are  quite  similar,  results  were  acceptable,  but  were  not 
the  best.  On  the  other  hand,  the  best  biorthogonal  wavelets 
corresponded  to  the  longest  filters  tested.  At  the  high 
compression  rates  that  were  used,  and  with  the  wavelets  that 
were  tested,  it  appears  that  the  number  of  vanishing  moments 
and  the  evaluation  of  regularity  are  of  secondary  importance 
when  compared  to  other  attributes  of  the  filter  banks. 

The  type  of  images  that  were  tested  was  selected 
purposely  in  order  to  test  the  wavelets  under  different 
situations.  In  the  first  instance,  the  axial  MR  images  were 
selected  in  order  to  test  the  compression  methods  under  the 
best  of  circumstances.  Since  most  of  the  image  information  is 
in  the  center  of  the  image,  and  the  periphery  is  mostly  dark, 
there  should  be  no  need  to  remove  border  artifacts,  and  the 
processing  algorithms  should  be  relatively  free  of  artifacts, 
and  indeed  this  is  the  case.  The  second  type  of  image 
corresponds  to  the  coronal  slices,  where  again,  most  of  the 
information  is  centered.  The  third  type  of  image  is  a  slice 
from  a  saggittal  study,  and  here,  since  there  is  a  larger 
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Figure  1.  Average  grades  for  the  four  images  compressed  with  the  wavelets  under  study. 
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amount  of  information  near  the  edges,  we  expected  results  except  on  the  saggittal  images  where  filter  b6.8  was  deemed 

with  more  reconstruction  artifacts.  The  final  image  is  a  to  be  the  best,  filter  b3.9  second  best,  and  b2.8  third  best, 

standard  256  x  256  Lena,  where  no  effort  was  made  to  avoid  Results  for  modal  analysis  are  very  similar.  At  a  33:1 

the  edge  artifacts.  compression  ratio,  wavelet  b2,8  is  graded  highest  in  all 

Averaged  results  show  that  for  the  case  of  a  33:1  instances  while  wavelet  3,9  received  the  same  grade  in  the 

compression  ratio,  grades  were  over  9,5/10  for  all  the  MR  case  of  axial  and  coronal  images.  Wavelet  b6,8  was  also 

modalities,  while  the  Lena  image  was  graded  with  a  9/10.  In  assigned  the  highest  grade  in  the  case  of  coronal  images.  In 

all  these  instances  the  wavelet  that  presented  the  best  results  the  case  of  a  50: 1  compression  ratio,  wavelet  b2,8  got  the 
was  b2.8,  except  for  the  case  of  saggittal  images,  where  highest  grades  in  all  cases  except  for  the  coronal  images, 

wavelet  b3.9  got  the  best  results.  For  a  50:1  compression  where  wavelets  b6,8  and  b3, 9  received  higher  grades, 

ratio,  results  show  the  same  pattern,  with  grades  over  8.4  for 
all  the  modalities,  with  filters  b2.8  giving  the  best  results 


Figure  2.  Results  for  a  50:1  compression  rate  and  the  best  filters  of  their  kind.  Top  row:  original  images; 
second  row:  biorthogonal  2,8;  third  row:  coiflet  C5;  last  row:  symlet  s8. 
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V.  Conclusion 

To  summarize  the  results,  we  can  say  that  among  the 
different  types  of  wavelets  that  were  tested  for  compression 
at  high  rates,  "symlets"  were  the  least  appropriate  for 
compression  of  images  in  general.  The  "coiflets"  that  were 
tested  in  this  study  were  similar  in  shape  to  most  other 
wavelets  in  the  study.  It  is  possible  that  the  degree  of 
assymmetry  in  the  filters,  while  small,  does  not  allow  an 
adequate  reconstruction.  In  the  case  of  the  "coiflets"  that  were 
tested,  the  results  were  almost  always  inferior  to  those 
achieved  by  biorthogonal  filter  banks.  In  some  instances  filter 
C5  outperformed  some  of  the  biorthogonal  filters,  but  was 
ranked  at  the  most  in  4th  place  among  the  filters  that  were 
tested. 

The  biorthogonal  filter  banks  were  the  best  performers  in 
the  group,  particularly  filters  b2.8,  b3.9  and  in  one  instance, 
filter  b6.8.  A  possible  explanation  for  this  performance  might 
lie  in  the  fact  that  these  are  the  longest  filters  among  those 
tested,  which  allows  a  better  representation  of  details.  These 
same  filters  obtain  the  best  results  for  the  qualitative 
evaluation,  as  well  as  for  several  quantitative  indexes  such  as 
the  Peak  Signal  to  Noise  Ratio  and  the  percentage  of  the 
energy  retained  after  compression,  and  are  similar  to  those 
described  in  [7]. 

It  appears  that  the  gross  wavelet  shape  as  described  by  the 
peak  to  peak  ratio  can  provide  a  general  idea  on  the 
performance  of  the  wavelets  but  is  not  a  reliable  index  to 
performance  in  image  coding  at  high  compression  rates. 
None  of  the  indexes  or  measures  that  have  been  used  can 
predict  individually  the  quality  of  results,  but  taken  as  a 
group  are  generally  consistent  in  the  evaluation  of  results.  It 
is  possible  that  a  combination  of  the  peak-to-peak  ratios  with 
the  percentage  of  the  energy  recovered  can  be  constructed  to 
form  an  index  that  might  predict  image  quality  after 
reconstruction  in  an  adequate  manner. 

Acknowledgment 

This  project  was  supported  by  CONACyT  grant  No. 
400200-5-29290. 


[5]  E.A.B  Da  Silva,  M.Ghanbari, 

"On  the  performance  of  linear  phase  wavelet  transforms  in 
low  bit-rate  image  coding,"  IEEE  Transactions  on  Image 
Processing,  vol.  5,  No.  5,  pp.  689-704,  1996. 

[6]  Abu-Rezq,  A.N.,  Tolba,  A.S.,  Khuwaja,  G.A.,  Foda,  S.G. 
"Best  parameters  selection  for  wavelet  packet-based 
compression  of  magnetic  resonance  images, "  Computers  and 
Biomedical  Research ,  vol.  21,  pp  449-469,  1999. 

[7]  Moreau  de  Saint-Martin,  F.,  Siohan,  P.,  Cohen,  A. 
"Biorthogonal  filterbanks  and  energy  preservation  property  in 
image  compression,"  IEEE  Trans.  Image  Processing,  vol.  8, 
no.2,  Feb.,  pp.  168-178,  1999 


References 

[1]  I.  Daubechies,  Ten  lectures  on  wavelets.  CBMS-NSF 
Regional  Conference  Series  in  Applied  Mathematics.  SIAM, 
Philadelphia,  PA,  1992. 

[2]  M.Antonini,  M.Barlaud,  P.  Mathieu,  I.  Daubechies, 

"Image  coding  using  wavelet  transform, "  IEEE  Trans.  Image 
Processing,  vol.  1,  pp.  205-220,  Apr.  1992. 

[3]  D. Field,  "Scale  invariance  and  self-similar  'wavelet' 
transforms:  An  analysis  of  natural  scenes  and  mammalian 
visual  systems,"  in  M.Farge,  J.C.R  Hunt,  and  J.C.  Vassilicos, 
eds.  Wavelets,  Fractals  and  Fourier  Transforms,  Clarendon 
Press,  Oxford,  1993. 

[4]  S.G.Mallat,  "Multifrequency  channel  decompositions  of 
images  and  wavelet  models,"  IEEE  Trans.  Acoust.,  Speech, 
Signal  Processing,  vol.  37,  pp  2091-2110,  Dec.  1989. 


Proceedings  -  23rd  Annual  Conference  -  IEEE/EMBS  Oct.25-28,  2001,  Istanbul,  TURKEY 


